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We have considered one dimensional Bose-Fermi mixture with equal densities and unequal masses 
using numerical DMRG. For the mass ratio of K-Rb mixture and attraction between bosons and 
fermions, we determined the phase diagram. For weak boson-boson interactions, there is a direct 
transition between two-component Luttinger liquid and collapsed phases as the boson-fermion at- 
■ traction is increased. For strong enough boson-boson interactions, we find an intermediate "paired" 

phase, which is a single-component Luttinger liquid of composite particles. We investigated correla- 
tion functions of such a "paired" phase, studied the stability of "paired" phase to density imbalance, 
and discussed various experimental techniques which can be used to detect it. 
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Recent developments in the cooling and trapping of atomic gases open exciting opportunities for experimental studies 
of interacting systems under well-controlled conditions. Using Feshbach resonances [l[ and/or optical lattices [1, Q it 
is possible to reach strongly interacting regimes, where correlations between atoms play a crucial role. The effect of 
interactions is most prominent for low-dimensional systems, and recent experimental realization [H H| of a strongly 
interacting Tonks- Girardeau (TG) gas of bosons opens new perspectives in experimental studies of strongly interacting 
systems in ID. In particular, one can experimentally study the behavior of Bose-Fermi (BF) mixtures Q in ID. Due to 
the lack of candidate systems in traditional solid state experiments, this topic did not attract much theoretical attention 
until recently. By now, properties of ID BF mixtures have been investigated using mean-field approximation [3], 
Luttinger liquid (LL) formalism [1, [^, [l(| [Hj], T- matrix approximation [T^], exact solutions and numerical 
methods [m EE EE The 

mean-field approximation is unreliable in ID, the LL approach describes the phase 
"-^J . diagram in terms of universal parameters which are hard to relate to experimentally controlled parameters, and the 
exact solution is restricted only to a certain region of parameter space. Most of the numerical work so far considered 
BF mixtures in optical lattices with fillings of the order of unity. In such a regime the analysis of the phase diagram is 
, ^ , ■ complicated, since the physics of the Mott transition plays an important role. In this article, we are mainly interested 
in the properties of ID mixtures without optical lattice, for the regime of parameters directly relevant to current 
£Nj ■ K-Rb experiments [l8j . Our main result is summarized in Fig. [TJ For K-Rb mixture with equal density of bosons 
and fermions, we find the "pairing" phase for moderate interaction strengths. Such a phase was previously discussed 
in Ref. Q for the case of almost equal masses and strong interactions based on LL formalism, but its relevance for 
^2 ' current experiments with K-Rb mixtures has not been addressed before. We note that the existence of "pairing" 
discussed in this article does not require the presence of commensurate optical lattice, as in Ref. [lj| . 
The general Hamiltonian of a ID BF mixture is given by 



where ty b , * / are boson and fermion operators, m b , mf are the masses, and 066, g b f are boson-boson and boson-fermion 
interaction strengths. Well away from confinement induced resonances [20j, ID interactions are given by 



i — i 

1 ofc j „ t .^b±^f±(m b + m f ) 

>T , 366 = 2nuj b± _ a bb and g bf = 2h ■ a b f , 

u b ±m b + uif±rrif 

where u> b ±,u>f± are transverse confinement frequencies, and a bb ,a b f are 3D scattering lengths. g bb and g b f can be 
controlled by changing transverse confining frequencies, or by varying scattering lengths using Feshbach resonances 
For a K-Rb mixture in the absence of magnetic field, a bb > and a b f < 0, so in this article we will study the regime 
when boson-fermion interaction is attractive and boson-boson interaction is repulsive. The phases of Hamiltonian ([TJ) 
in the most general case depend on four dimcnsionless parameters, which we choose to be m b /mf, n b /rif, 

m b g bb ^/m b m f g b} 

766 = —2 and 76/ = 2 (2) 

n n b n ^JnynJ 

Here rif,n b are fermion and boson densities, and 766 and 76/ are dimensionless interaction parameters. Similar to Lieb- 
Liniger model (2lT |. strongly interacting regime corresponds to 766, 17&/I 3> 1- If 766 < (attractive bosons), the system 
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FIG. 1: Phase diagram for K-Rb mixture with equal density of bosons and fermions. 7^5 and 76/ are defined by Eq. ((2j), and 
we show here the attractive side (7^ < 0). Different symbols and colors identify the three different phases involved: squares 
are in the two-component Luttinger liquid (LL) state, diamonds are in the collapsed state, whereas circles stand for the paired 
phase we found in between. Sample points, considered in Figs. [2] [3] and |4j are identified by stars. Solid line is the mean-field 
prediction Q, which should be good for small jbb and matches the data for smallest ybb available. Dashed lines are tentative 
boundaries between the different phases. Notice the existence of a threshold in jbb to get the paired phase. 



is always unstable towards boson collapse. For "f bb > the system can still collapse for 7^ < 0, or phase separate for 
76/ > (H . If these two scenarios are not realized and densities of bosons and fermions are incommensurate, then 
from LL theory [Id [13 

one expects a two-component LL, with power law decay of all correlations. If densities are 
commensurate, one can expect [g, [ll[ a nontrivial pairing, resulting in the exponential decay of certain correlation 
functions and in the opening of the gap. In what follows we will concentrate on the latter case for an experimentally 
relevant K-Rb mixture, so we will fix n b /nf — 1, m b /mf = 87/40 and consider 7^/ negative. 

For numerical purposes we consider the discretized version of Hamiltonian ([T]) to be an open boundary chain with 
unity lattice constant and L sites. Similar to the Lieb-Liniger model being a low filling fraction limit of Bose-Hubbard 
model [22j, Hamiltonian (JTJ is the low filling fraction limit of the following lattice Hamiltonian: 



Hi 



L-l L 
i=l i=l 



t b (b\b i+l +H.c)-t f (flf i+x +H.c)) +J2(^btb i (blb i -1) + U bf btb i flfi) ■ (3) 



We note that many different lattice Hamiltonians give continuum Hamiltonian ([T]) in the low density limit, and the 
choice of lattice Hamiltonian is not unique. For low fillings v b and Vf, dimensionless interaction parameters are given 
by 

~ Ubb , Ubf ,.s 

^ ~ 2^ and ™ * 2^m^- (4) 

Most of our simulations were performed at densities v ~ 1/4, but some of the results were checked for v ~ 1/8. The 
fact that we use finite filling fractions only slightly changes Eq. (j4|) , but does not affect the phase diagram qualitatively. 
We stress here that our results were checked also against commensurability effects by testing them at filling fractions 
v = 23/96 and v = 25/96 too. Thus we can safely say that pairing effects are not depending on the particular v value 
chosen. The expectation values (. .J) of one- and two-body operators over the ground state of TIl have been evaluated 
by means of the DMRG method [23| , which provides a practically exact solution for any value of the couplings and 
allows one to measure correlation functions with both statistics on the equal footing. We used up to L — 96 chains, 
with local dimension d = 10 (up to 4 bosons per site) and truncation up to m = 256 states. Discarded probabilities 
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amount to less than e = 5 • 10~ 7 . Small values of jbb with weak interparticle interactions, however, are not easy to 
study with this method, since high occupation number for bosons should be taken into account. We thus resorted to 
mean-field predictions in this area. 

Local density profiles were measured and are plotted in Fig. [5] for two sample points (corresponding to stars in 
Fig.[l}. In our simulations we calculated the following correlation functions: bosonic Green function Gf,(i,j) = (b\bj), 
fermionic Green function Gf(i,j) — (fjfj) and "pairing" correlation function Gbf(i,j) = (AT A 3 -) = (b\fjfjbj), 
where the "pairing" operator is described by Aj = fjbj. Their values have been plotted in Fig. [3] for the same sample 
points as before. Fourier transforms of bosonic and fermionic Green functions give respective momentum distributions, 
while Fourier transform of Gbf{i,j) is related to momentum distribution of composite particles. Such momentum 
distributions are shown in Fig. [4j In order to calculate them, it is crucial to choose the properly defined free-particle 
eigenmodes due to open boundaries 4>k(J) = v^/C^+l) si n (fc j) with k — nj^(n = 1, ...,L). In addition to these 
correlation functions, density-density correlations D a> p = (n a (i)np(j)) — (n a (i)) (np(j)) were measured as well. Here 
a and j3 can take any value from {b, /}. For finite size simulations, subtraction of the non-connected part is necessary, 
since open boundary conditions give density profiles with Friedel oscillations. Examining of density correlations 
without this substraction could lead to misinterpretation of results. 

To extract the long-range behavior of correlation functions, we restricted the analysis to the regions far away from 
the boundaries. In order to check whether a certain correlation function has a power-law or exponential decay, it is 
sufficient to test the simplest power-law and exponential forms 

G*(x)= £ ,. am j sin(w*2; + (^), (5) 

where the oscillating term is absent for the pure bosonic Green function. The fermionic oscillation frequency is correctly 
given by the fermionic density ujf — nvf ~ tt/4. The exact value is given by the density in the system bulk, which is 
slightly larger due to open boundary conditions. For convenience, let us introduce the Luttinger parameters Kb, Kf 
and Kbf, which are related to a^, o.j and o.bf in Eq.(f5]) as o.b — l/(2Kb), otf(bf) = (Kf(bf) +l/-^/(6/))/2- Equation ([5]) 
gives the asymptotic form of the correlation functions in the thermodynamic limit far from the boundaries. To 
quantitatively extract the Luttinger parameters K* for finite L, one has to take into account carefully the effects 
of open boundary conditions (OBCs). We refer the reader to the detailed analysis of Ref. [HI and recall here that 
G*(i,j) will depend on chord functions di,(x) = d {x\2{L + 1)) = (2(L + sin (ttx/2(L + 1)) of all 2i, 2j,i—j, i+j. 

The form of the correlation function which needs to be fitted to extract Luttinger parameters are given by (Gbf is 
modified in the same way as G/) 

G b (i,j) cx [d L (2i)d L (2j)]*k [d L (i + j)d L (i-j)]-^k , (6) 
G f (i,j) <x [d L (2 i )d L (2j)f iiK, -^7 ) Ml(*+jK(*-j)H (K/+7 ^ (7) 
{A, [d L (i + j)] K ' (-lf^ri sin(7rw(t-i) + tp Q ) + A x [d L (i-j)] K f S m(iru(i+j) 

We start our tour around the phase diagram by looking at small attractions between bosons and fermions. The 
sample point jbb — 6.0 and 7h/ = —1.36 is considered in the first panels of Figs. [2] and [H Looking at density profiles in 
Fig. [21 one can notice that both the atomic species spread out around the whole lattice, exhibiting Friedel oscillations 
due to hard walls at the boundaries. Wings are cut off for the sake of plot clearness. Due to attraction between 
bosons and fermions such oscillations are in-phase, but the values of the two densities differ on the order of the second 
digit. Coming to correlation functions (illustrated in Fig. [3J first panel), one can easily recognize a power-law decay 
for both the bosonic Green function Gb(i,j) (black circles) and the (oscillating) fermionic Green function Gf(i,j) 
(red squares). "Pairing" correlation function Gbf(i,j) exhibits power-law behavior as well. Thus this phase is a two- 
component Luttinger liquid. Such phase has two gapless sound modes, and all correlation functions have algebraic 
decay. 

Increasing the interspecies attraction for small "/bb will lead to a collapse. More precisely, bosons form a small 
region with high density where fermions will be attracted up to Pauli-allowed density v = 1. Existence of such 
maximal density is an artifact of our lattice discretization, and is not expected in the absence of a lattice. According 
to mean- field theory [7j, the first order phase transition between the two component LL and the collapsed phase 
should take place as boson-fermion attraction is increased for any value of the boson-boson interaction. Within mean- 
field theory transition line is given by 7^ = jbb^ 2 , and it is shown in Fig. [1] as a solid line. The result of mean 
field calculation agrees well with the data set for smallest 7^ considered. For large jbb, mean- field calculation is not 
expected to give an adequate description of the system, and for sufficiently large "fbb and attractive fermion-boson term 
76/ < the system belongs to a third intermediate phase, see Fig. [1] for a sketch. The population distributions and 
the correlations for a sample point in this intermediate region of parameters are plotted in the second panels of Figs. [2] 
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FIG. 2: Density profiles: black circles are bosons, red squares fermions. In the first panel we show results for a typical point 
in the usual mixed phase (755 = 6.0,76/ = —1.36). Both the species are spread around the chain and no pinning is evident 
between densities. In the second panel, the "locking" of densities for the paired phase is evident. Right panel corresponds to 
the values 766 = 6.0,76/ = —8.14. Friedel oscillation frequency is in both cases given, as expected, by u» = nv* with v* ~ 1/4 
being the species' density in the bulk of the system. 



and [31 A strong locking of one density profile on top of the other is the most striking feature in Fig. [D^b). Indeed, not 
only the Friedel oscillations are in phase as they were in Fig. [2a), but the difference between boson and fermion local 
densities is bounded to be less than 10" 4 in the bulk, which is two orders of magnitude smaller than in the case of 
two-component LL. Furthermore, the strong locking of the two densities suggest that a composite particle made by a 
boson and a fermion, polaron, could be the new elementary object to look at. In the second panel of Fig. [3] all three 
types of correlation functions are plotted: G b (black circles), G f (red squares) and Gbf (green diamonds). In contrast 
to two-component LL, single species Green functions clearly exhibit an exponential decay with a correlation length 
of few sites. However, the "pairing" correlations Gbf still decay slowly as a power-law. Taking the open boundary 
conditions into account as described before, we get a Luttinger parameter Kbf — 0.95 ±0.02. As shown in Fig.[U such 
a dramatic change in the decay properties of correlation functions is witnessed by momentum distributions of the two 
species (and the composite one). Indeed, the Fermi step of individual fermionic atoms is no more there as it is in the 
case of two-component LL, and also the once peaked Bose distribution is considerably spread out now. In contrast, 
we highlight that momentum distribution of paired composite particles clearly exhibit a Fermi step around kbf = it /4 
consistently with the filling. Thus, this phase can be understood as the "paired" phase of bosons and fermions. The 
existence of such a "paired" phase has been predicted in Ref. Q based on LL theory arguments for mixture with 
equal masses for large 7^. Indications of the existence of such phases have also been briefly presented in Ref. [13], 
but the phase diagram has not been studied in detail. Figure [T] presents the phase diagram for the K-Rb mixture, 
and shows that "paired" phase can be realized for moderately strong Bose-Bose interactions. Boundaries between 
different phases were determined comparing algebraic and exponential fits of single species correlation functions, and 
by observing the "locking" of one density profile on the top of the other. 

Looking at density-density correlations, one can address another distinctive feature of "paired" phase. As predicted 
in Ref. Q, oscillating part of all three density correlation functions D a ^{x) (b — b, f — /, b — f) decay with distance 
with the same algebraic exponent 

D*(x)\ 27r v ~ \x\~ r sm(ujx + tp). (8) 

The frequency of oscillations is twice the particle density uj ~ Itxv ~ n/2. As pointed out by Ref. Q, exponent r 
should be intimately related to the Luttinger parameter for paired particles Kbf, i.e. r = 2Kbf (we note that in our 
notations Kbf = K+/2, where K + in defined in Ref. Q). Thus parameter Kbf can be extracted independently from 
D a ,/3, using the fitting procedure which takes OBC into account. We checked that all density-density correlations 
decay with the same exponent, and extracted value of Kbf equals 0.97 ± 0.02. This is in good agreement with the 
pairing correlation fits (see before) which give the value 0.95 ± 0.02. Both these derivations have been checked for 
scaling with respect to both the system size L (from 48 to 128) and the DMRG truncation parameter m (up to 
320, discarded down to < 10 -8 ). The values of Kbf has been thus confirmed to survive to the thermodynamical 
limit. Based on all evidence, we can unambiguously state that we have shown the existence of the "paired" phase 
predicted by Cazalilla and Ho @, even with moderate interactions and unequal masses of the two atomic species 
for the experimentally relevant case of K-Rb mixture. Such a phase survives up to the thermodynamical limit with 
an almost unchanged exponent near 1. We recall here that the universal prediction of unity value holds only at the 
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transition point itself which, on the other hand, is difficult to be precisely addressed due to its BKT nature. 




FIG. 3: Correlation functions: same sample points as in Fig. [21 same color code. For clarity, oscillations are not shown 
and only the decay of the envelope functions is presented. G(x) means here G(L/2, L/2 + x). Green diamonds stand for 
composite particles created by A 1 " = b' f' . In the first panel (two-component Luttinger liquid state) all the three types of 
correlation functions exhibit an algebraic decay. For bosons Gb(x) oc a; _1// ' 2Ki> ' with Kb = 1.45 ± 0.05, whereas for fermions 
Gf(x) oc a; _:L /( 2if /) _if // 2 with Kf = 0.98 ± 0.02. On the other hand, for "paired" phase (second panel) only Gbf shows an 
algebraic decay ~ £ _1 /( 2 - ftr f>/)--R"b// 2 w ;th the Luttinger parameter Kbf = 0.95 ± 0.02, while Gb and G/ decay exponentially 
with distance. 




FIG. 4: Momentum distributions: same sample points and same color code of correlations as in Fig. [3] In the first panel, data 
for the two-component Luttinger liquid phase, where a Fermi surface for individual fermions and a tightly peaked distribution 
for bosons can be clearly seen. In the second panel, momentum distributions in the "paired" phase. We highlight the washing 
out of the Fermi surface for individual 40 K and the wide broadening of the bosonic 87 Rb distribution. On the contrary, a sharp 
step- like feature in the composite particles' mode occupation is present around 7r/4, indicating the algebraic decay of "pairing" 
correlation function Gbf(x). Bosonic and fermionic particle correlations decay exponentially (see the inset, which is the same 
data as in the second panel of Fig. |3j but in log-linear scale). 

The stability of "paired" phase against the population imbalance between the two species was also briefly studied 
(see Fig. [5]). It turns out that for small enough density differences, the "locking" of the densities survives in the sides 
of the box, whereas in the middle a peak or a hole arises in the bosonic profile. The system phase separates into 
the region which exhibits "pairing" and the region with unequal densities. For larger density imbalance, the system 
becomes unstable to collapse, and the "paired" phase is washed out. This example qualitatively illustrates that the 
"paired" phase can be observed even if the densities of bosons and fermions are not exactly equal, but the imbalance 
is smaller than some threshold. More detailed studies of phases with unequal densities lie beyond the scope of this 
work. 

Let us now concentrate on possible experimental techniques to detect the "paired" phase. One notable feature, 
which distinguishes the "paired" phase from the two-component LL phase, is the presence of the gap for out-of-phasc 
density modes. Deep in the "paired" phase, the energy scale for the gap is set by Fermi energy ~ Tr 2 h 2 n 2 /(2m). 
The presence and the size of the gap can be measured using RF spectroscopy [25|, [26j]. Energies of sound modes 
can be measured using Bragg scattering [I?!, and the "paired" phase has only one sound mode, as opposed to 
two-component LL phase, which has two modes. Disappearance of the out-of-phase sound mode also affects qual- 
itatively the response of the system to the variations of the external potential, since out-of-phase collective mode 
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FIG. 5: Density profiles in the presence of an extra boson (first panel) or fermion (second one) with respect to filling 1/4 in 
the paired region (same settings as Fig. [2| . The two species' profiles retain their "pairing" in the wings and exhibit a peak or 
a hole in the bosonic density in the center of the box. For larger density imbalance, the system becomes unstable to collapse. 



in external trap is much higher in energy compared to in-phase mode. Fourier transforms of bosonic and fermionic 
correlation functions, (nj,(k)) = J 4>k{x)Gb(x)dx and (rif(k)) — J <pk(x)G f(x)dx, are given by momentum distri- 
butions, shown in Fig. 3] They can be measured using ballistic time-of-fiight experiments, since during ballis- 
tic expansion momentum distributions get mapped into real space densities [Bj]. The momentum distribution of 
composite particles (iibf(k)) — J <fik(x)( , &l(x)' ( &'f(x)' , &f(0) 1 $>b{0)}dx has a strong Fermi step, and can be written as 
(ribf(k)) = J (nb(k — p)nf(p))dp/(2n). We note that it is different from J (rib (k — p))(rif (p))dp/ (2tt) , thus the presence 
of the "pairing" results in nontrivial noise correlations in time-of-fiight images [28( . Finally, we would also like to point 
out the method to measure the correlation functions based on interference of two independent ID clouds [29|, [30(. For 
bosons (fermions) the average of the square of interference signal (\A b ^)(L)\ 2 ) of two segments of length L is related 

to an integral of the Green's function as (|A b (j)(L)| 2 ) = J Q L dxdyG b (f)(x — y) 2 . The interference signal appears at 
wave vectors Qb(f)> which depend on masses of interfering particles. If one measures (\Ab(f)(L)\ 2 ) as a function of L, 

then in principle dependence of Gb(j)(x) on distance x can be extracted, since Gb(f)(L) 2 — \ — b Q^i — ■ The same 
technique can be used also to measure Gbf{x), but in this case the information will be contained in the oscillations 
of the product of Bose and Fermi densities at wave vector Qb + Qf. Since in "paired" phase Gbf(x) decays much 
slower with distance than Gb(x)G f(x), "paired" phase will be characterized by strong correlations in the fluctuations 
of bosonic and fermionic interference fringes. 

To summarize, we have considered one-dimensional Bose-Fermi mixture with equal densities and unequal masses 
using DMRG. For the mass ratio of K-Rb mixture and attraction between bosons and fermions, we determined the 
phase diagram, which is shown in Fig. [T] For weak boson-boson interactions, there is a direct transition between 
two-component Luttinger liquid and collapsed phases as the boson-fermion attraction is increased. For strong enough 
boson-boson interactions, we find an intermediate "paired" phase, which is a single-component Luttinger liquid of 
composite particles. We investigated correlation functions of such a "paired" phase, studied its stability to density 
imbalance, and discussed various experimental techniques which can be used to detect it. 

We thank E. Demler, V. Gritsev, R. Fazio, and F. Dolcini for useful discussions. This work has been developed 
using the DMRG code released within the "Powder with Power" project (www.qti.sns.it) 
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